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ABSTRACT 

We compute the probability distribution for the spin of a black hole following 
a series of minor mergers with isotropically distributed, non-spinning, inspiraling 
compact objects. By solving the Fokker-Planck equation governing this stochas- 
tic process, we obtain accurate analytical fits for the evolution of the mean and 
standard deviation of the spin distribution in several parameter regimes. We com- 
plement these analytical fits with numerical Monte-Carlo simulations in situations 
when the Fokker-Planck analysis is not applicable. We find that a ~ 150 M Q 
intermediate-mass black hole that gained half of its mass through minor mergers 
with neutron stars will have dimensionless spin parameter \ = a/M ~ 0.2±0.08. 
We estimate the effect of the spin of the central black hole on the detection range 
for intermediate-mass-ratio inspiral (IMRI) detections by Advanced LIGO and 
extreme-mass-ratio inspiral (EMRI) detections by LISA. We find that for real- 
istic black hole spins, the inclination-averaged Advanced-LIGO IMRI detection 
range may be increased by up to 10% relative to the range for IMRIs into non- 
spinning intermediate-mass black holes. For LISA, we find that the detection 
range for EMRIs into 10 5 M & massive black holes (MBHs) is not significantly 
affected by MBH spin, the range for EMRIs into 10 6 M & MBHs is affected at the 
< 10% level, and EMRIs into maximally spinning 10 7 M MBHs are detectable 
to a distance ~ 25 times greater than EMRIs into non-spinning black holes. The 
resulting bias in favor of detecting EMRIs into rapidly spinning MBHs will play 
a role when extracting the MBH spin distribution from EMRI statistics. 

Subject headings: black hole physics — gravitational waves 



1. Introduction 



A growing body of evidence from observations, numerical simulations, and comparisons 
between the two, suggests the existence of a population of intermediate-mass black holes with 
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masses in the M ~ 10 2 — 1O 4 M range (e.g., (IMiller fc Colbertll2004j ) and references therein) 



These intermediate-mass black holes may capture compact obj e cts (stellar-mass black ho 



or neutrons stars) and me rge with them ( Taniguchi et al. 



2000; 



Mouri fc Taniguchil l2002al ibl; iGiiltekin. Miller, fc Hamilton 



Miller fc Hamilton 



2004 



2006 



2002a 



O'Leary et al. 



OS 



b; 



2006 



Mandel et al.l 120071 ). In addition to adding to the black- hole mass, the merging compact 
objects will also contribute their orbital angular momentum to the spin angular momentum 
of the central black hole, leading to the evolution of the black-hole spin through a sequence 
of such minor mergers. 

We might expect the typical spin of a black hole to be low if a significant fraction of 
its mass has been added via minor mergers with compact objects whose angular momentum 
at plunge is distributed isotropically. The angular momentum imparted to the black hole of 
mass M by a compact object of mass m is L ^ oc mM. (We include only the orbital angular 
momentum, not the spin angular momentum of the compact object, since the latter is lower 
than the former by a factor of order m/M, which we assume to be small for minor mergers.) 
This causes the dimensionless spin parameter of the hole x = S\/M 2 = a/M to change by 
~ L bj/M 2 oc m/M. After iV ~ M/m such mergers, necessary for the hole to grow to mass 
M, the typical dimensionless spin parameter of the hole will be x oc {m/M)yfN ~ yJm/M. 



As discussed by IMiller! (120021 ) and iHughes fc Blandfordl (12004 ). the angular momenta 
of black holes that grow through minor mergers undergo a damped random walk. The 
damping comes about because retrograde orbits, which subtract angular momentum from a 
black hole, plunge from a last stable orbit (LSO) at a higher radius than prograde orbits, so 
more angular momentum is subtracted following retrograde inspirals than is added following 
prograde ones. 

In this paper, we make an analytical approximation to the spin change induced by a 
minor merger and s olve the Fokker-Planck equa tion to obtain the evolution of the spin prob- 
ability distribution (IHughes fc Blandfordl 120041) . (We use a sim pler one- dimensional version 
of the Fokker-Planck equation than lHughes fc Blandfordl (120041 ). since we are interested only 
in the evolution of the magnitude of the spin, not its direction.) We find that for black holes 
with x ^ ym/M , the spin x evolves proportionally to M~ 2m as the mass grows via minor 
mergers (rather than M~ 2 , which would be the case without damping). We determine the 
asymptotic values of the expected mean of the spin distribution and its standard deviation 
in the limit of infinitely many minor mergers: x ~^ 

■sJl.hm/M and a -> y/0.7m/M. We 
also describe the evolution of the spin distribution in other parameter regimes, e.g., when 
y/m/M > x > m/M. 

Our Fokker-Planck analysis fails when the mass ratio m/M is not sufficiently low, so 
for those cases we resort to Monte-Carlo numerical simulations. We find that if the mass 
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of the central black hole grows from M = 5m to M = 10m by capturing five objects of 
equal mass m, the mean spin of the resulting black hole is \ ~ 0-5, nearly independent of its 



initial spin (iMillerl (120021 ) obtained similar results). However, if the central black hole grows 
from M = 50m to M — 100m (e.g., a M = 70 M black hole growing to M — 140 M by 
capturing fifty m = 1.4 M & neutron stars), its resulting spin is rather low, \ ~ 0.2 ± 0.08. 

The combination of the spin of the central black hole and the inclination of the in- 
spiraling object's orbit can have a significant effect on the gravitational-wave signal from a 
low-mass-ratio inspiral. We compute the increase in the Advanced-LIGO detection range 
for intermediate-mass-ratio inspirals (IMRIs) due to the spin of the central black hole. We 
find that the detection range, averaged over orbital inclinations, may increase by ~ 3 — 10% 
relative to the range for inspirals into non-spinning black holes for the expected values of 
black hole mass and spin. We provide an approximate expression for the dependence of the 
Advanced-LIGO IMRI detection range on spin [see Eq. fl24")) ]. We also compute the change in 
the LISA extreme-mass-ratio-inspiral (EMRI) detection range due to the spin of the massive 
black hole. We find that the range for inspirals into M = 10 5 M black holes is nearly 
independent of their spin, because the frequency at the last stable orbit (LSO) is away from 
the minimum of the LISA noise curve. On the other hand, the inclination-averaged detection 
range for IMRIs into rapidly spinning M = 10 7 M black holes is ~ 25 times greater than 
into non-spinning ones. The detection volumes are proportional to the cube of the range. 
This will create a bias in favor of detecting inspirals into rapidly spinning black holes, which 
in turn will have consequences for the extraction of massive-black-hole spin function from 
LISA EMRI statistics. 

The paper is organized as follows. In Sec. [2, we provide the background for our calcu- 
lation of the spin evolution via minor mergers. In Sec. [3], we describe analytical solutions 
of the Fokker-Planck equation for spin evolution. In Sec. HJ we describe Markov-Chain nu- 
merical simulations of spin evolution. In Sec. 0, we evaluate the dependence of the detection 
ranges for low-mass-ratio inspirals averaged over orbital inclination angles on the spin of the 
massive body, in the context of both Advanced LIGO and LISA. 



2. Spin evolution 

We assume that the distribution of the orbital inclination angle i relative to the central 
black hole's spin is isotropic at capture. Here i is defined via 

cost= — z , (1) 
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L z is the object's orbital angular momentum in the direction of the black hole's spin, and Q 
is the Carter constant. We further assume that the inclination angle i remains approximately 
constant over the inspiral (jHughesI l2000l ) . so the distribution of inclinations at the LSO is 
also isotropic, Pr(cost) = 1/2. 

In the low-mass-ratio limit, the amount of angular momentum radiated in gravitational 
waves during the plunge and ringdown is smaller by a factor of ~ m/M than the angular 
momentum at the LSO. Therefore, we assume that the merging object contributes its orbital 
angular momentum at the LSO to the angular momentum of the black hole. The spin of the 
black hole after a minor merger, x', is related to the original spin % y i a 



X 



V( X M 2 + L,) 2 + Q, 



(2) 



(M + m) 2 

where m is the mass of the small object, M is the mass of the hole, and we assume m -C M. 

The constants of motion L z and Q at the LSO can be obtained as a function of t by 
demanding that the potential R and its first and secon d derivatives in r are zero at the LSO 
(see Chapter 33 of (IMisner. Thorne. fc Wheelerlll973l )): 



R 
R 



[E(r 2 + x 2 M 2 ) - L zX M] 2 - (r 2 - 2Mr + x 2 M 2 ) [mV + (L z - xME) 2 + Q] 



dR 

dr 



d 2 R 
dr 2 



at LSO. 



(3) 



It is possible to make analytic approximations to the values of L z and Q at the LSO 
based on appropriately averaging the analytically known con stants of motion at the LS O 
for prograde and retrograde equatorial orbits (cf. Eq. (9) of (IHughes fc Blandfordl 120041 ) ). 
In particular, for \ <C 1, the plunging object's dimensionless "total angular momentum" is 
given by 



Mm 



Mm VI 2 



where we correct a mistake in Eq. (4) of jMilleril2002 h. Then L z and Q follow from Eq. ([!]): 



2 \ 3 



3/2 



X COS L 



(4) 



L z = cos i a/ L\ + Q; Q = sin L 2 + Q. 



(5) 



3. Fokker-Planck equation for spin evolution 



The black-hole spin evolution is a stochastic process. The probability distribution func- 
tion of a stochastic process, however, can be described the deterministic Fokker-Planck 
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equation: 



0_ 

Of 



d_ 

dx 



\p>(x,t)f(x,t)} + 



i a 2 



[a 2 (x,t)f(x,t)] , 



(6) 



2dx 2 



where /i = (dx)/dt is the mean drift and a 2 = ((dx) 2 )/dt is the stochastic variance. In this 
Section, we derive approximate analytical solutions to the Fokker-Planck equation in several 
interesting parameter regimes. 

For simplicity, assume that all merging objects have the same mass m. We parametrize 
the mass of the black hole by a dimensionless "time" parameter t = M/m. The change in 
the spin x after a merger follows from Eq. ([21) : 



We can compute L at plunge as a function of x an d cost by solving Eqs. ([3]), then 
substituting the result into Eq. ([7]) to obtain dx as a function of t, x, and cost. Although 
this process is simple in principle, such a numerical computation makes it impossible to 
obtain analytic expressions for (dx) and ((dx) 2 ), which are necessary if we wish to solve the 
Fokker-Planck equation. (Here, brackets denote averaging over cost.) 

We could, of course, try to obtain empirical analytic fits to the numerical solutions for 
(dx) and ((dx) 2 ), but it turns out that there is a simpler approach. The approximate formula 
for L given in Eq. (J3J) is valid only when ^ < 1; when % ~ 1, Eq. (j3J) overestimates L by 
as much as 40%. Remarkably, however, using this incorrect approximation for L in Eq. (J7J) 
generally yields very accurate expressions for (dx) for a wide range of x- So long as xt 3> 1 
(i.e., x ^ ui/M), an expansion of Eq. ([7]) to the first order in l/(xt) yields the following 
simple analytic expression for the mean drift in x : 



This expression is accurate to about 1% for all values of x so l° n g as ^ 10. Similarly, the 
analytic expression for the stochastic variance of the spin is 






8 




(9) 



This expression underestimates the variance by > 10% for very high spins, but is generally 
accurate to a few percent for lower spins which are expected as a consequence of minor 
mergers in the Advanced LIGO setting. 
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We can now substitute Eqs. (jHJ) and 
bility evolution ([6]) to obtain 



d_ 
dt 



f(x,t) 



into the Fokker-Planck equation for the proba- 



(10) 




This is a one-dimensional equation unlike the three-dimensional equation derived in (IHughes fc Blandford 



2004 ). since we choose to focus on the evolution of the magnitude of the spin, not its direc- 



tion. Still, this is a rather complicated equation that does not easily separate. Fortunately, 
for many applications it is not necessary to solve the complete equation. 

Equation (jTOl) was derived under the assumption ^ > 1. If we further assume that 
X 2 t 3> 1 (i.e., x ^ a/ m/M, then the mean spin evolution is dominated by 



dx 



X 

dt t 



where a = -2 - 4^/9 « -2.63. (This result can also be obtained directly from Eq. (jHJ).) 
Thus, the mean spin evolves according to 



X 



Ao 



Mo 
M 



2.63 



(12) 



(compare with Eq. (26) of ( IHughes fc Blandfordll2004l ). where the exponent is approximated 
by 2.4). 

If the assumption x 2 t 3> 1 is not satisfied, and instead x 2 t "C 1, but \t 1 so that 
Eq. ffTU]) still holds, the evolution of the probability function may be approximated as 



df(t,x) 
dt 



_o_ nf(t, x ) 

dx V Xt 2 



+ 



1 d 2 



( W, x) 



2d X 2 V t2 



(13) 



This equation can be solved by separation of variables: f(t,x) — T(t)X(x), where the 
solution for T is T(t) = exp(—k/i), X is the solution to 



2 X 2 X"-A X X + AX-kx 2 X 
and A; is a constant. The mean spin grows roughly as 



0. 



X 



2 

T 



2 
t' 



(14) 



(15) 
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so after t > 2to (i.e., after the black hole captures half its mass via minor mergers), x 2 t ^ 1- 

The spin growth and spin decay terms in Eq. (fTUl) cancel when the spin is approximately 
equal to 

nr /l5 , \ 

^V^ W VT (16) 



(Compare with iMillerl ( 20021 ). who estimated the mean spin to be \j2^jJmjM) = \/2/t 
based on numerical simulations.) 

We can estimate the second moment of the proba bility distribu tion by approximating 
the solution to Eq. ( 1101) by a Gaussian (as suggested by Miller J2003)): 



f(x,t) 



2na 



exp 



(x-mr 

2a 2 (t) 



(17) 



(A Gaussian turns out to be a good approximation except at small x, when the tails at 
X > X are larger than those at x < X-) Substituting this Gaussian into Eq. (fTDl) . keeping 
only the lowest-order terms in tx, and setting x — Xi we obtain 



1 da 
o dt 



a 
1 



t 2 a' 



il + bx 2 ), 



where b = 4^/2/9 — 1. If a 2 t ^> 1, then a evolves in the same way as x when % 2 t ^> 1: 

2.63 



(J W cr 



O"0 



Mo 
M 



(19) 



What if cr 2 t < 1? This might be the case of interest if, say, the initial spin of a black 
hole created during some process is known precisely, and we wish to estimate future spin 
evolution through minor mergers. In this case, the second term on the right-hand side of 
Eq. (fTS|) dominates, and if x is small or does not change significantly, o grows according to 



a 



id • i>\ 2 ) ( r - 7 J 



(20) 



In either case, a asymptotes to the solution 



a 



2(1 + b X 
—at 



(21) 



For large t, a ~ ^2/ (-at) w ^0.7/*; Miller! (I2002f ) estimated tr to be ^(m/M)/V2 
sj\j (2t) based on numerical simulations. 



- 8- 



Lastly, consider the case when < 1. In this case the orbital angular momentum of 
the plunging object is comparable to the spin angular momentum of the black hole, and 
Eq. (fTUj) is incorrect, since it was derived under the assumption xP ^> 1. If the black hole is 
initially non-spinning or has spin x ^5 however, a single minor merger will bring its spin 
to x ~ yVl/t according to Eq. (J7J). This case can be treated with a Monte-Carlo numerical 
simulation as described in the next section. 



4. Spin evolution via Monte Carlo simulations 

We have carried out Monte Carlo simulations of spin evolution through minor mergers 
in order to confirm the analytical estimates presented above, based on the approximate 
Fokker-Planck equation. Our simulations also allow us to access the small-t regime where 
the Fokker-Planck approach is not valid, but where our physical approximations for low- 
mass-ratio inspirals still hold. Since these simulations were performed numerically, there 
was no need to make analytical approximations to d\ following a merger; instead, we solved 
Eqs. (J3J) directly and obtained dx via Eq. ([7]). 

In Figure [1] we plot the spin distribution of a black hole of mass t = M / m — 10 that 
started out with either spin x = 0-1 or X = 0-9 at i = M/m = 5 before growing via minor 
mergers. This corresponds, for example, to an intermediate-mass black hole that grows from 
M = 50 M Q to M = 100 M & by capturing m = 10 M & black holes. The distributions for 
both values of initial spin are roughly Gaussian, although with shorter-than-Gaussian tails 
(we plot the actual Monte-Carlo histogram for the \ = 0.9 case for comparison with a fitted 
Gaussian). We see that for these small values of t, the initial value of the spin is largely 
forgotten after the black hole captures half of its mass through minor mergers. The means 
of the spin at t — 10 are x = 0-49 for the initially slowly-spinning hole and x = 0-51 for the 
initially rapidly-spinning hole. The standard deviations at t — 10 are a = 0.17 for initial spin 
X = 0.1 and o = 0.18 for initial spin x — 0.9 (the initial standard deviations are zero in both 
cases, i.e., the initial spins are presumed to be precisely determined). These results agree 



with Fig. 1 of (jMillerl 120021 ). Because the values of t involved are so small, the Fokker-Planck 
equation (flUj) does not apply: at t = 5, the angular momentum of the inspiraling object at 
the LSO is comparable to or larger than the spin angular momentum of the black hole even 
for large initial black hole spins. 

In Figure [2] we plot the spin distribution for a black hole of mass t = M/m = 100 
that started out at t — M/m = 50 at either spin x — 0.1 or x — 0.9 before growing via 
minor mergers. This corresponds, for example, to an intermediate-mass black hole that 
grows from 70 M & to 140 M by capturing M = 1.4 M neutron stars. The means of the 
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Fig. 1. — Monte-Carlo predictions for the black-hole spin distribution following black- hole 
growth via minor mergers from t = M/m = 5 to t = M/m = 10. The histogram shows the 
spin distribution at t = 10 for a black hole with initial spin x — 0.9, and the solid curve is a 
Gaussian fit to that distribution. The dashed curve is a Gaussian fit to the spin distribution 
at t — 10 for a black hole that has initial spin x — 0.1 at t — 5. 
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spin at t — 100 are x = 0.162 for the initially slowly-spinning hole and x = 0.233 for the 
initially rapidly-spinning hole. The final spin in the initially rapidly-spinning case decreases 
as X ~ Xo(V^o) -2 ; rather than x ~ Xo(^Ao) -2 ' 63 as predicted by Eq. f[T2l . That is because 
the spin begins to approach the asymptotic value of x ~ \/l.5/t ~ 0.12 as predicted by 
Eq. (|T6l) . and the rate of spin evolution decreases because x 2 t is no longer much greater than 
one. The initially slowly-spinning case does not quite satisfy %t 3> 1, so the Fokker-Planck 
analysis is suspect; however, Eq. (|T5|) . relevant since x 2 t < 1 in this case, provides a roughly 
accurate estimate of spin growth. The standard deviations at t = 100 are a = 0.066 for 
initial spin x — 0.1 and a = 0.084 for initial spin x — 0.9; the predicted asymptotic value 
of the standard deviation according to Eq. (j2ip is cr = 0.087. The mass ratios considered in 
this paragraph may be plausible for intermediate-mass-ratio inspirals into interm ediate-mass 



black holes that would be detectable with Advanced LIGO (IMandel et al.l 120071 ). 



Finally, we perform a Monte-Carlo simulation of the evolution of a spin distribution from 
t = 1100 to t = 1200 where the starting mean spin is x = 0.72 and the starting standard 
deviation is a = 0.016. In this case, x 2 t ^ 1 holds throughout the evolution, so this example 
can be viewed as a test of our Fokker-Planck analysis. Based on Eq. (TI21) . we expect the 
spin at t — 1200 to decrease to x — 0.57; in fact, we find x(t = 1200) = 0.58. Since a 2 t <C 1, 
we expect the standard deviation to grow via Eq. (j20|) to a = 0.022 at t — 1200; in fact, 
a{t = 1200) = 0.021. 

The Fokker-Planck analysis should give excellent results in the regime of very large t, 
such as those corresponding to minor mergers of stellar- mass compact objects with ~ 10 6 M Q 
massive black holes in galactic centers. (The extreme-mass-r atio inspirals preceding su ch 



minor mergers are an interesting class of potential LISA sources ( lAmaro-Seoane et al.ll2007l ).) 
On the other hand, if a large range of t must be covered, Monte-Carlo simulations become 
expensive. Thus, the Monte-Carlo numerical methods and Fokker-Planck analysis can be 
viewed as complementary techniques. 



5. Effect of black-hole spin on detection ranges for low-mass-ratio inspirals 

The frequency of the last stable orbit before plunge is strongly influenced by the black- 
hole spin and the orbital inclination. Prograde inspirals into rapidly spinning black holes 
will have much higher LSO frequencies than inspirals into non-spinning black holes or polar 
inspirals into spinning black holes of the same mass, while retrograde inspirals into rapidly 
spinning black holes will have lower LSO frequencies. For example, for a maximally spinning 
Kerr black hole, the frequency of the LSO of a retrograde equatorial inspiral is twice lower 
than for a polar orbit, while the LSO frequency of a prograde equatorial inspiral is six 
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Fig. 2. — Monte-Carlo black-hole spin distribution following black hole growth via minor 
mergers from t = M/m = 50 to t = M/m = 100. The spin distribution for a black hole with 
initial spin x — 0.9 is shown with a solid curve, and one for initial spin x = 0.1 is shown 
with a dashed curve. 
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times higher than for a polar orbit. Even for a more moderately spinning black hole with 
X = 0.4, there is almost a factor of two difference between LSO frequencies for prograde and 
retrograde inspirals. 

The signal-to-noise ratio (SNR) for the detection of gravitational waves from inspirals 
depends on where the LSO frequency falls on the noise power spectral density curve of the 
detector. Although some inclination angles will increase SNR and others will decrease it, 
we might generally expect that average detection range for inspirals into spinning black 
holes will be higher than into non-spinning ones. ("Average" refers to averaging over the 
isotropically distributed orbital inclination angles of the inspiraling object.) This is because 
of the cubic dependence of the detection volume on detection range, which is proportional 
to SNR: if, say, 10% of all inspirals have their SNR boosted by a factor of three, these will 
be seen three times further and the detection volume for these kinds of inspirals will go up 
by a factor of 27, so the average volume in which detections can be made will increase by a 
factor of ~ 3, and the average detection range will grow by the cube root of 3. 

Conversely, this average detection range increase can manifest itself as a bias in favor 
of detecting inspirals into rapidly spinning black holes rather than slowly spinning ones. 
Thus, a numerical estimate of the detection range increase due to black hole spin is useful 
for determining whether a high fraction of rapidly spinning black holes among detected 
inspirals is an indication of the prevalence of such black holes in the universe, or whether 
this is merely a selection effect. 

We use the simple scaling 

\kf) 2 \ oc r 7/ " (22) 

for the frequency-domain gravitational wave. The square of the signal-to-noise ratio p 2 is 
proportional to 

p2 K /" /max \kf) 2 \ df K /" /max r^ df (23) 

Jfmm Sn(f) J / min S n (f) 

Here, S n (f) is the noise power spectral density of the detector, / max is the frequency of 
gravitational waves from the last stable orbit, and f m i n is the low-frequency cutoff for the 
detector for Advanced LIGO, where / min = 10 Hz, or the frequency of gravitational waves 
one year before plunge for LISA. We set / max equal to twice the orbital frequency at the 
LSO, which we obtain numerically as a function of the black-hole mass M and spin x an d 
of the orbital inclination angle cost by solving Eq. (j3J). 

The distance to which an event can be seen is proportional to SNR, p, so the detection 
volume is proportional to p 3 . Therefore, we average p 3 , computed via Eq. (j2"3"l) . over the 
different inclinations cost (uniformly distributed through the range [—1, 1]) in order to 
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compute the expected increase in the detection volume for a given values of Xi an d then take 
the cube root to compute the increase in the average detection range. 

We have computed detection ranges for A dvanced LIGO using this method with the 



noise power spectral density S n (\f\) taken from (jFritschelll2003l ). Fig. [3]shows our computed 
ratio between (i) the average Advanced-LIGO detection range for intermediate-mass-ratio 
inspirals into black holes of a given mass and spin and (ii) the detection range for IMRIs into 
Schwarzschild black holes with the same mass. For low spins x ^ 0-4, which are typical for 
intermediate-mass black holes of ~ 100 — 200 solar masses that gained a significant fraction 
of their mass via minor mergers, we can approximate the detection range increase due to the 
inclusion of central black hole spin as 

R-w* .. 1 + tf("). (24) 



Range no _ spin \100 M Q / 

This is the ratio of detection ranges; the ratio of detection volumes is a cube of this ratio. 

The effects of cosmological redshift are not significant for Advanced-LIGO IMRIs when 
the black-hole spin is small. Even prograde equatorial inspirals of neutron stars into M = 
100 M Q black holes spinning at x — 0-9 are only detectable to z m 0.2 at an SNR threshold 
of 8. The cosmological redshift has the same effect as increasing the black-hole mass, so 
including redshift increases the ratio of detection volumes at higher spins. For the purposes 
of including redshift in Fig. [3], the inspiraling object mass was set to m — 1.4 M and a 
detection threshold of SNR = 8 was assumed. 

The results described here do not include higher-order (m ^ 2) harmonics of the orbital 
frequency. Higher harmonics are not significant when black-hole spins are small, since in 
that case they affect both the spinning and the non-spinning rates roughly equally, and so 
the ratio does not change. However, for high values of spin, the ratios would probably drop 
somewhat relative to those given in Fig. [31 since including higher-frequency harmonics would 
contribute more to increasing the detection range for inspirals into non-s pinning holes than 



into rapidly holes with prograde orbits (cf. Fig. 6 of (IMandel et al.l 120071 )). 



We also compute the dependence of the LISA EMRI detection range on the massive 
black hole spin. We consider EMRIs of m = 10 M objects into M = 10 5 M , M = 10 6 M , 
and M = 10 7 M massive black holes. We assume that a detection is possible at an SNR 
threshold of 30. (Setting the threshold to 15 changes the results at the 10 — 20% level.) 
Cosmological redshift must be included for LISA EMRIs since they can be seen to z ~ 1 — 2. 
This means we must specify the inspiraling object mass and the SNR detection threshold, 
since these are necessary to determine the cosmological redshift of the most distant detectable 
source. 
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Fig. 3. — The ratio between the inclination-averaged Advanced-LIGO detection range for 
intermediate-mass-ratio inspirals into Kerr black holes of a given spin and the detection 
range for IMRIs into non-spinning black holes. The solid curve represents black holes with 
mass M = 100 M ; the dashed curve, mass M = 200 M Q . 
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LISA EMRIs only sweep through a fraction of the frequency band during the observation 
time. Therefore, f m \ n for LISA is set not by the detector threshold, but by the frequency of 
the gravitational waves emitted one year before plunge. We compute f min by evolving the 
gravitational-wave frequ ency back in time from plunge for one year using the prescription of 



Barack fc Cutler! (J200J) (Eqs. (28) and (29)). 



For M = 10 5 M Q , the spin of the black hole is almost irrelevant: once we average over 
orbital inclinations, the spin affects the detection range at a level of at most a few percent. 
This is because at these low masses, most of the SNR comes from the portion of the inspiral 
at much higher radii than the LSO, so the exact freque ncy of the LSO does not p lay a very 
significant role (cf. Fig. 8 and associated discussion in (lAmaro-Seoane et al.l 120071 )). 

Figure H] shows the dependence of the average EMRI detection range on the massive- 
black- hole spin for M = 10 6 M & . The average detection range for EMRIs into rapidly 
spinning black holes of mass M = 10 6 M & is ~ 13% larger than for EMRIs into non-spinning 
black holes. For M = 10 7 M , the detection range for EMRIs into rapidly spinning black 
holes is increased by a factor of ~ 25 over those into non-spinning black holes, as shown in 
Fig. This greater sensitivity to black hole spin is expected, since for these massive black 
holes most of the SNR comes from the cycles near the LSO. However, this should not be taken 
to mean that inspirals into rapidly spinning M = 10 7 M black holes are likely to dominate 
LISA EMRI observations. Figures H] and [5] show detection range ratios only; the inclination- 
averaged detection range for an EMRI into a maximally spinning M = 10 7 M black hole is 
actually less than the detection range for an EMRI into a non-spinning M = 10 6 M & black 
hole. On the other hand, this large ratio does mean that there is a strong detection bias 
in favor of rapidly spinning black holes, which must be taken into account when statistics 
of EMRI observations are inverted to gather information about the massive-black-hole spin 
distribution. 
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Fig. 4. — The ratio between LISA detection ranges (at SNR= 30) for extreme-mass-ratio 
inspirals of m — 10 M compact objects into Kerr black holes of mass M = 10 6 M and a 
given spin vs. non-spinning black holes. 
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Fig. 5. — The ratio between LISA detection ranges (at SNR= 30) for extreme-mass-ratio 
inspirals of m — 10 M compact objects into Kerr black holes of mass M = 10 7 M and a 
given spin vs. non-spinning black holes. 
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